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Abstract 

The immunity of microwave-induced magneto-resistance oscillations and corresponding zero re- 
sistance regions to the direction of (circular) polarization of microwave is studied in this paper. 
We propose that a spontaneous circular motion of the whole electron fluid would stabilize the 
system and minimize the polarization sensitivity of the oscillatory DC resistance. Results of a 
self-consistent calculation capture the qualitative features of the experimental observation. 
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The observation of "Zero-(DC)-Resistance State" (ZRS) in high mobility two-dimensional 
electron systems (2DES) under magnetic field and microwave (MW) radiation presented a 
surprise to the physics community [l[ • When the system, under crossed uniform magnetic 
field and small dc bias, is irradiated by microwave of sufficient intensity, the longitudinal 
(dissipative) resistance R xx develops strong oscillatory dependence on the magnetic field. At 
low temperature and high radiation power the minima of the oscillations evolve into ZRS. 

Intense interests in the physics community have been aroused by this unusual non- 
equilibrium phenomenon. Theoretical efforts aimed at identifying the microscopic mech- 
anisms accounting for the MW-induced resistance oscillations (MIRO) and ZRS. Early the- 
ories employed the picture of photon-excited transport assisted by short-range scattering 3, 
Q, 01, and later the scenario of MW-induced oscillations in the nonequilibrium electron dis- 

n 

tribution function was argued to be the leading cause of MIRO[6[]. Other mechanisms based 
on different physical pictures were also proposed^, 0, 0, Q. As for ZRS, though there exist 
proposals that do not invoke negative resistance (NR) instability 11], the pattern formation 
model [3], which is based on NR, remains the most popular. On the experimental side, 
the outstanding issues are activated temperature dependence with large energy gaps [131 ]. 



immunity of MIRO and ZRS to the polarity of circular polarization of MW 



14J , suppression 



of MIRO and ZRS by in-plane magnetic fields 151 ]. and multi-photon processes in the ZRS 



formation 



161 ]; most of these findings are not readily accommodated by existing theories. 



Among other issues, the polarization immunity of DC resistance poses a particularly 
challenging test. Experiment by Smet et a/.jjjl established that MIRO are insensitive to 
the polarization state of the MW radiation. However, transmission data also shows that 
active cyclotron resonance absorption occurs only when the circular polarization matches 
the magnetic field orientation (denoted as CRA), although ZRS is observed under both 
MW polarities. This observation was inconsistent with the two most prevailing theories, the 



IB, 3 



and the 



impurity and/or phonon assisted inter-Landau level (LL) transitions model 
non-equilibrium distribution function scenario [6J. These theories predict oscillations with 
correct period and phase, however, with substantially different amplitudes for the different 
polarizations; the MW photoconductivity of the cyclotron resonance inactive (CRI) state 
is smaller by a factor ((a; — u c ) 2 + T 2 ) / ((u + u c ) 2 + T 2 ), where T" 1 is a phenomenological 
lifetime. The factor reflects the huge difference in the (AC) Drude conductivity for opposite 
circular polarities. 
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Polarization immunity of the DC resistance posts a big challenge to the understanding of 
the ZRS. To produce the "same" MIRO and ZRS, there must be an additional mechanism 
that compensates for the discrepancy in the energy absorption rate between the CRA and 
CRI states. In this paper, we propose that a spontaneous circular motion of the whole 
electron fluid previously proposed by Ng and Dai 17j could stabilize the system and minimize 
the polarization dependence of the oscillatory photoconductivity. The spontaneous circular 
motion arises in the cyclotron resonance (CR) favorable orientation whenever the intensity 
of incident MW radiation exceeds a threshold value. 

To see how this could happen, we examine the general transport equation for the center 
of mass (CM) coordinate R(t) of the electron liquid where the effect of impurity is included 
to the second order 

t 

m R(t) = (-e) \ g(t) + ^R(t) x Bj + aV 5(t) J dt'x (R(t) - R{t'); * - f ) , (1) 

— oo 

where a = rii\u\ 2 /n, rii is the density of impurity, |ii| 2 indicates the strength of the impurity 
potential, n = N/V is the carrier density, and x is the (equilibrium) retard density-density 
response function without the MW term. The equation is derived in the CM frame where 
the MW field is eliminated and the electron liquid sees moving impurities following the path 
~ —R(t). After impurity averaging, we obtain Eq. 17). We note that similar equation 
has also been proposed by Lei et a/.j^]. The electric field relevant to this phenomenon is 



En. 



of form E(t) = E\ cos (uit) + E% sin (fit) + Ed, where E2 = ±z x E\ and 
The first two terms represent a circularly-polarized MW radiation with frequency uj and E& 
stands for a small DC bias. The plus/minus sign in the expression of E2 indicates CRA/CRI 
state. For small DC bias R(t) can be written as R(t) = RAc(t) + Roc{t), where Rac{P) is 
the (dominant) part induced by MW field and Rr>c(t) ~ vt is a small DC correction. 

Eq. flTJ suggests that the linear DC resistance is determined completely by Rac{P) since 
there is nothing else in the equation. As sketched in Fig.l (left solid circles), RAc(t) for 
opposite polarization states have very different amplitudes of motion in the Drude model, 
and to explain the polarization immunity, this big difference in R,Ac(t) must be minimized. 
A plausible way for this to occur is to have a spontaneous circular motion in the CRA 
direction for both CRA/CRI states. If this spontaneous circular motion dominates over the 
Drude motion, the difference in Rac(1) between the CRA and CRI states will be minimized. 
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In reality the composite trajectory of the spontaneous and Drude motion can be rather 
complicated (Fig.l right-hand solid traces) and a self-consistent numerical calculation has 
to be performed to test this idea. 



CRA 



CRI 




O 



FIG. 1: Schematic trajectory of the AC motion of the electron fluid in the absence (left) and 
presence (right) of the spontaneous circular motion (dashed circle) for CRA (upper) and CRI 
(lower) states. 

To test our proposal, we study a trial solution of Eq. (JTJ with 

R(t) = v-t+R f (t)+R s (t) = v-t+Ai cos (ut + 5)+Bi sin [pt + 5)+di cos {u z t)+0i sin {u) t t) , 

(2) 

where Rf(t) describes the CM motion directly coupled to the radiation field and R s (t) 
depicts the spontaneously generated oscillatory mode with frequency We shall call them 
"fast" and "slow" modes in the following in view of their frequency difference which will be 
shown numerically later. 5 is the phase delay of the "fast mode" with respect to the AC 
driving. Notice that we keep only the base harmonic modes in R/(t) and R s (t). This can be 
justified for the "fast" mode since in the weak radiation limit the size of the circular orbit 



R c ~ y A\ + B\ will be much less than the magnetic length l c = ^Jc/eB (We set h = 1 in 
the following) which is quite long for the B field (~ IT) relevant to this phenomenon 171] . 
Notice that such an argument does not apply to the "slow" mode whose amplitude is not 
governed by the radiation strength. 

— # 

We next expand the impurity-scattering induced damping force F^t) = 
aVa { j J dt'x (R(t) — R(t');t — t'j in a Bessel- Fourier series and keep only the base- 

— oo 

frequency oscillating terms. By putting the approximate R(t) and Fj(t) into Eq. (1) and 
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compare, we obtain three force-balance equations 



= -eE d - - v x B + F v , (3) 

c 



mR f (t) = -eE{t) - - R f {t) x B + (4) 
mR 8 (t) = -~R s (t)xB + Fj s \t). (5) 



Here 

» » OO 

F v = -a d d qq ^ Im X (Q, Q- v + mu + nu^ J 2 t (z (q)) J 2 t (y (q)) (6) 

m,n=— oo 

is the time-averaged damping force, where -2(g) 2 = (g • Ai) 2 + (g ■ Bi) 2 and y(q) 2 = (q • 
«i) 2 + (q • Pi) 2 - For the "fast mode", F^ f) (t) ~ Hi cos (tut + 5) + G 1 sm(iut + 5) where 
Hi = 7r (1) ■ Ai - 7r (2) ■ Bi and G% = 7r (1) ■ B x + 7r (2) ■ A 1 whereas for the "slow mode" , 
F^(t) ~ hi cos + sin (w r t), where /ii = -n^' -ai — n^' ■ j3i and #1 = n^' • /3i+ir^' -di. 
The 7r's are given by 



^(1) a 
n 2 

a 



— / / d d g g 2 Rex (<f, <f • u + ttuo + nco> 7 

' m,n=— 00 

» » OO 

7r^ = — / / d d g g 2 /mx (g, g • u + mcj + nujj 



j m (z(m j; " m) 



m 

m,n=—oo 



z(q) 
J m {z{q)) 



J 2 n(y(o)l 



z{q) 



JMo)), (7) 



and 7rW, tc^' can be obtained from tt^\ tt^ by simply interchanging z <-> ?/, m <-> n and 
a) <-> We shall consider the limit w — > when calculating 7r's in the following since we 
are interested in the linear-response DC current. 

After some algebra we obtain the equation of motion for "fast mode" 



muo 



-a 



(1 T 00J00) + vrWf + [tt^] 2 • A x = eE , (8) 



where — /+ corresponds to CRA/CRI state. The equation is of Drude form representing an 
oscillation driven by AC electric field under impurity-induced friction tt^ that depends on 
the amplitudes of both "fast" and "slow" modes (z(q) and y(q))- A corresponding reactive 
correction 7r^ also appears which can be interpreted as a mass correction [17]. The phase 
delay 6 is given by tan 5 = / [ma) 2 (i =pw c /w) Notice the explicit polarization 

dependence of A\ in Eq. (8). 

Since there is no external driving force, the spontaneously-generated "slow" mode is 
determined by the self-sustainability requirement that the corresponding frictional force 
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vanishes, i.e. 



7T 



(2)' 



0. (9) 



Putting it into Eq. (5) we obtain another equation for "slow mode" , 

— mw ; 2 (1 — uJc/uj;) a± = tt^ ■ ol\. (10) 

The two equations determine self-consistently the amplitude and frequency of the "slow 
mode" . Notice that the spontaneous circular motion is always in the CR-favorable direction, 
independent of the polarization of the MW radiation. Polarization dependence enters only 
indirectly through z(q) which appears in both tt^' and n^' . 

Eq.s (8)-(10) form a set of self-consistent equations determining the "fast" and "slow" 
mode amplitudes A\, ot\ and the "slow" mode frequency u) z . These equations are solved 



numerically. The longitudinal DC resistance R xx = — lim 

i? — s-0 



v~ 2 E 



is computed afterward. 



We have employed the density-density response function of non-interacting electron gas 



in constant magnetic field [18j in our calculation, with Landau levels (LLs) broadened phe- 
nomenologically into Lorentzians, i.e. S(e — nu c ) — > n^T/ ((e — nuj c ) 2 + V 2 ). Since exper- 
imentally Co ~ 200GHz is fixed while sweeping the B field, we use uo as the basic unit and 
set T ~ Cu, r(T) ~ 0.2a), E F ~ 10a; and keep 20 LLs in our calculation, consistent with 
the low field (B ~ IT) and intermediate temperature (T ~ IK) setting in experiment. We 
shall vary magnetic field B and use the frequency ratio u n = u/uj c as abscissa in present- 
ing our results. We also define normalized radiation intensity I N = (^eEo/m*u! 2 l) 2 where 
I = y/ uj c /uj l c (m* = 0.068m is the effective mass of conduction band electrons in GaAs) as 
well as normalized amplitudes c = A\j I, d = a%/ 1 and renormalized frequency x = io I juj c 
in our calculation. 

To further simplify calculation and analysis, we make another approximation of keeping 
only \n\, \m\ < 1 terms (zero- and one-photon processes) in the impurity induced forces in 
our calculation; this is consistent with keeping only base harmonics in our trial trajectory. 
Correspondingly, we restrict our calculation to the frequency range u n = uj/uj c ~ 1-2 since 
MIRO occur at the weak B field side Q > uj c and higher frequency range will be dominated 
by transitions between higher LLs (\n\, \m\ > 1 processes). The approximation of keeping 
only single photon process in the "fast mode" can be justified in the weak radiation field 
limit. However the approximation can not be justified a prior for the spontaneous (slow) 
mode. We have estimated the effect of multi-photon processes associated with the "slow 



mode" on our self-consistent equations and found that it does not affect terms associated 
with Imx(q, q-v + ncUj+Cu) but would enhance terms associated with Imx(q, q-v + nujj) and 
Imx(q, q • v + uuj i —uj). To mimic these effects we introduce a correction factor a > 1 that 
multiplies the latter two terms in n^' . The same factor is introduced to the counterparts 
in 7rW and F v (R xx ) to ensure consistency in our calculation. As long as a does not deviate 
from a = 1 too much, the qualitative behavior of the solutions is not sensitive to a. We shall 
present our calculation results with a = 1.2 in the following. 

Numerically we find that we may roughly divide the region u n ~ 1-2 into negative 
resistance (NR) region that centers about u> n ~ 1.2, and positive resistance (PR) region that 
peaks at u n ~ 1.7. For radiation intensity (I N ) t h < I N < (I N )muiti-u>, where (7^)^^0.001 
numerically (corresponding to E ~ 7 V/cm) is a threshold value above which "slow" mode 
appears and (I N ) mu iti-u> ~ 0.02 (E Q ~ 32 V/cm) is a value where multi-photon processes of 
MW radiation become important, we find that the solutions for c, d and x fall within the 
ranges listed in Table 1. 



polarity 


frequency range 


d 


c 


X 


<CRA> 


NR(w„ : 1.1-1.5) 


0.1-0.2 


0.2-0.4 


~ 0.3 




PR(w n : 1.5-1.9) 


0.2-0.3 


0.2-0.4 


~ 1.1 


<CRI> 


K : 1-1-2.0) 


~ 0.1 


0.03-0.08 


~ 0.3 



TABLE I: The numerical range of the self-consistently determined d, c and x for CRA and CRI 
states in the regime uj n ~ 1-2. 

Notice that the amplitude of "fast mode" (c) is much smaller in the CRI state when compared 
with CRA state as expected. However, the amplitudes d of the "slow mode" in the two 
polarization states are comparable. The low-frequency "slow mode" (cUj ~ 0.3u; c ) appears 
in the whole frequency range of CRI state and the NR region of CRA state. Unexpectedly, 
d 7^ solutions also exist at the PR region of CRA state with rather large amplitude and 
near-CR frequency. 

The effect of the "slow mode" on the polarization-dependence of DC resistance and ZRS 
can be roughly understood as follows: when the incident radiation is sufficiently strong and 
NR instability shows up, a spontaneous circular motion is generated in the electron fluid 
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in CR favorable direction. Mathematically, low-frequency fictitious photons represented 
by n 7^ processes are spontaneously generated and additional photon-assisted transport 
channels open up. As a result normal dissipation is enhanced and NR is cured. This 
mechanism dominates the whole frequency range of CRI state and the NR region of CRA 
state. The spontaneous circular motion also turns the original fast-rotating CRA orbit 
(left-hand solid circle of Fig.l) into the "lace" of a slowly-rotating orbit (right-hand 
dashed circle). The strong mixing of the two orbits destroys the identity of well-defined 
Landau orbit and suppresses photon-assisted scattering in the CRA state, thus drawing the 
two states more close to each other. Mathematically, the suppression comes through the 
increase in the argument (y(q)) of Bessel functions. 




FIG. 2: DC resistance versus frequency ratio for CRA and CRI states under two values of radiation 
intensity: (a) I N = 0.005 (moderate intensity) and (b) I N = 0.01 (high intensity). (R) — with slow 
mode included in calculation, (r) — slow mode is excluded. 

We show in Fig. 2 the calculated DC resistance versus frequency ratio both with (R) and 
without (r) "slow mode" included in our model. Notice the different plotting scale for the two 
calculated resistances in the figure. We see that when the radiation is of moderate intensity 
I N = 0.005(i?o ~ 16 V/cm) (Fig. 2a) , NR region becomes pronounced in the absence of the 
"slow mode" for both CRA and CRI states and the disparity in the oscillatory amplitude 
of the two polarization states is large. 

The entering of the "slow mode" suppresses the oscillatory amplitudes, especially for CRA 
state, which effectively reduces the amplitude discrepancy between the two polarizations; 
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meanwhile the NR instability is almost completely healed by the u^-photon-excited inter- 
LL transitions in our highly simplified model. As the radiation strength further increases 
to I N = 0.01 (£0 ~ 22V/cm) (Fig.2b), oscillation becomes stronger for CRI state but is 
approaching saturation for CRA state in the absence of the "slow mode". In this case, the 
suppression effect plus the photon-assisted transport processes associated with the "slow 
mode" make the oscillatory shape of the two polarization states even closer and lift the NR 
region to ZR for both polarities. We note also that at the B field region near CR where 
significant absorption takes place, no d 7^ solution is found in our calculation and the 
oscillation amplitude of CRA and CRI curves are considerably different, in agreement with 
experiment 14 ] . 

In conclusion, we propose and demonstrate numerically that a spontaneous slowly- 
rotating circulating current in a 2DES under magnetic field will be generated when the 
incident circularly-polarized MW radiation exceeds certain threshold intensity. The sponta- 
neous mode can cure the NR problem and provides a plausible explanation for the observed 
ZRS and polarization immunity of the DC resistance. Our calculation is crude because of 
the many approximations we made and the results are only in semi-quantitative agreement 
with experiment. Nevertheless, we believe our theory has provided a promising starting 
point to understand the physics behind the ZRS phenomenon and associated polarization 
immunity. 
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